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ABSTRACT 

Frequency noise in the variations of the Earth’s obliquity (tilt) can modulate the 
insolation signal for climate change. Including this frequency noise effect on the incoming solar 
radiation, we have applied an energy balance climate model to calculate the climate fluctuations 
for the past one million years. Model simulation results are in good agreement with the 
geologically observed paleoclimate data. We conclude that orbital noise in the Earth system may 
be the major cause of the climate fluctuation cycles. 

INTRODUCTION 

The cause of the 100-Kyr ice age cycles during the past one million years has been a 
scientific mystery despite decades of studies. Although the Earth’s orbital eccentricity signal 
varies with a 100-Kyr cycle, the resulting changes in solar heating are too small to be 
climatically significant. How can such a strong climate response arise from such a seemingly 
weak forcing? Numerous hypotheses and models in terms of internal feedbacks such as ocean 
circulations and atmospheric C0 2 concentrations have been proposed to understand this mystery. 
However, the origin of the climate fluctuations remains a perplexing and enduring problem. The 
purpose of this paper is to demonstrate that the cause of the major climate fluctuation cycles is 
mainly due to the frequency noise in the Earth’s obliquity variations, which is distinct from the 
commonly considered signals in the orbital eccentricity and internal feedbacks. 

CLIMATE RESPONSE TO ORBITAL NOISE 

Climate response to insolation signal is intrinsically nonlinear in a noise environment of 
the Earth’s orbital motion, including variations in eccentricity, obliquity and precession. The 
insolation amplitude signal alone, which is due to orbital variations, is unable to reproduce major 
spectral peaks as seen in the observed climate records. However, noise from the frequency 
variations in the obliquity may amplify the small signal of insolation forcing and cause the 
climate to fluctuate from glacial to interglacial conditions. The frequency noise from the 
obliquity is generated by the solar-luni and planetary gravitational torques exerted on the 
equatorial bulge of the Earth [1-6]. 

NOISE-SIGNAL COOPERATIVE PHENOMENON 

Based on the convolution theory of the non-stationary time series of orbital motion, the 
frequency noise in the obliquity can be expressed by [4,5]. 
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where the amplitude A„ frequency y, and phase for i = 1,2,3.... 80 are astronomical constants 
[7]. Liu [2-5] and Liu and Chao [6] have shown that the signal of the incoming solar radiation 
has been embedded in the frequency noise N(t) of the obliquity. This insolation signal S(t) is 
governed by [8]. 

S(t) = AR s Ae(t) + mA(esino)) (2) 

where A R s is the difference of the insolation received at the top of the atmosphere for 1 ° change 
in obliquity e(t) e sin W is the precession variation and m is a function of the duration of the 
tropical year, the orbital eccentricity and the latitude. The cooperative phenomenon between the 
signal S(t) and noise N(t) is shown in Fig la for the past one million years in which the two 
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Figure 1. Cooperative phenomenon between insolation signal S (t) and frequency noise N (t) of 
the Earth’s obliquity variation. 
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dashed threshold lines for S(t) represent maximum and minimum values of insolation deviation at 
latitude 65°N. Hence, amplitudes of S(t) are subthreshold [9]. To seal eN(t) as a pulse function, 
the pulse magnitudes (length of the bars in Fig lb) are chosen such that they can be normalized 
between ± 2 with zero mean. Fig lb shows a bipolar pulse modulation train of the insolation. 

Each time the magnitude of S(t) plus its corresponding N(t) crosses the threshold a pulse can be 
written to the time series of the insolation. It has been found that the timing of the insolation 
pulses coincides with the occurrence of major paleoclimate changes [4,5]. The pulse modulation 
train for the insolation signal S(t) as shown in Fig lb can be expressed by [10] 

c(t)=Za.P('- t .)-Sa j P(t-t i ) (3) 

I j 

where a 7 and a y are the magnitude of the insolation pulses and f, and tj are time indices of rapid 
melting and prolonged glaciation. 


CLIMATE FLUCTUATION 


The noise-signal cooperative phenomenon may cause pulses in the insolation signal S(t) . 
Including this pulse modulation effect, equation of the O-D energy balance climate model 
becomes [3-5] T " 
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Figure 2. Temperature fluctuation simulation T(t) from model equation (4) and the oxygen 
isotop 5 I8 0 concentration variations from [12]. 
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where d, d 2 d } T 0 , A and B are geophysical constants [11]. We have performed a series of 
computations to test model response to the pulse-modulated insolation. The fluctuation of the 
Earth’s surface temperature T(t) in equation (4) for latitude 65°N’is shown in Fig. 2, in good 
agreement with variations in the oxygen isotop content 8 ,a 0 of deep-sea sedimentary cores as 
obtained by [12], Extension of model simulation reveals that the current warming trend of the 
climate is almost over and will give way to forecast for a small ice age. 

DISCUSSSION AND CONCLUSION 

A faithful representation of climate fluctuations in computer form is needed to clarify 
events of the past, better to predict the future. But the route to that end still lies ahead of us 
because the cause of the major climate fluctuation cycles remains a scientific mystery today 
despite decades of investigation. This paper has revolutionized the knowledge of orbital 
dynamics of the Earth system and climatology. It has been found that climate response to 
insolation signal is intrinsically nonlinear in a noisy environment of orbital motion. The 
insolation amplitude signal alone, which is due to orbital variations, is unable to reproduce major 
spectral peaks as seen in the proxy climate data. The frequency noise in the Earth’s obliquity, 
however, may produce a noise-signal cooperative phenomenon that amplifies the small 
insolation forcing. This new forcing mechanism could enhance the probability of jumping 
between the glacial and interglacial climate states. 
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